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ABSTRACT 

Analytical investigation of dynamic stall on HAWT 
rotor loads was conducted. Dynamic stall was modeled 
using the Gormont approach on the MOD-2 rotor, treat- 
ing the blade as a rigid body teetering about a fixed 
axis. Blade flapwise bending moments at station 370 
were determined with and without dynamic stall for 
spatial variations in local wind speed due to wind 
shear and yaw. The predicted mean flapwise bending 
moments were found to be in good agreement with test 
results. Results obtained with and without dynamic 
stall showed no significant difference for the mean 
flapwise bending moment. Ihe cyclic bending moments 
calculated with and without dynamic stall effects were 
substantially the same. None of the calculated cyclic 
loads reached the level of the cyclic loads measured 
on the MOD-2 using the Boeing five-minute-average 
technique. 

NOMENCLATURE 

a Speed of sound 

Lift coefficient 
c Blade chord 

I Teeter moment of inertia 

M Flapwise bending moment 

M Jfech number 

r Distance along blade 

R Blade radius 

s Distance measured normal to axis of rotation 
v Blade of deflection 

V Wind speed 

w Aerodynamic load per unit span 

W Local relative velocity 

a Time rate of change of angle of attack 
a Angle of attack 

Dynamic stall parameter 
6 Yaw angle 

§3 Teeter axis tilt 

Teeter angle 

y Blade mass per unit span 

9 ^faster blade position angle 

£2 Rotor angular velocity 

INTRODUCTION 

Large horizontal axis wind turbine rotors experience 
on the order of ICr stress cycles during a projected 
30-year lifetime. Cyclic loads due to gravity, tower 
shadow, wind shear, and yaw contribute to stress 
cycling as well as wind turbulence which can be viewed 
as a cyclic disturbance in the reference system of the 
blade. With the exception of gravity, the loads 

mentioned above are aerodynamic in nature and, since 
dynamic stall has been Identified as playing a signif- 
icant role in influencing the performance and fatigue 
life of the vertical axis wind turbine [l], it has 
been conjectured that dynamic stall plays a similar 
role in influencing cyclic loads of horizontal axis 
wind turbines. 

The primary characteristics of dynamic stall are the 
occurrence of stall beyond the angle of attack associ- 
ated with the static stall, the dependence of the 


loads on the time rate of change of the angle of 
attack as well as the angle of attack, and hysteresis 
in the lift, drag, and pitching moment. These charac- 
teristics would be expected to Increase both mean and 
cyclic loadings. Noll and Ham [2] in their investiga- 
tion of the effects of dynamic stall on small wind 
systems concluded that dynamic stall could increase 
fixed pitch horizontal axis wind turbine normal loads 
and moments by "about 10 percent." Accordingly, this 
study was undertaken to examine the possible effects 
of dynamic stall on the MOD-2 wind turbine. Extensive 
data on mean and cyclic flapwise loads have been 
obtained for the MOD-2 [3] . 

ANALYSIS 

Geometry 

The analytical model used to investigate the effect of 
dynamic stall on the mean and cyclic flapwise bending 
moments consisted of a rigid two-bladed rotor teeter- 
ing about a fixed point. An existing program [4] was 
modified to treat a teetering rotor. Figure 1 shows a 
blade element, the wind, Vy, and the coordinate 
system. Looking upwind, the rotation of the rotor is 
counterclockwise with the rotor teeter angle 
designated as i[>. The axis of rotation, x, is fixed in 
space so that rotor yaw, represented by the angle 5 , 
Is due to shifting of the wind rather than motion of 
the yaw axis. Each blade of the two-bladed rotor lies 
in the same plane so that there Is no precone to the 
rotor. 

The blade positions differ by an angle of 180° and the 
teeter angles of the two blades differ in algebraic 
sign. One blade is designated as the master blade and 
the position angle, 6 , and the teeter angle, 4 >, are 
given for the master blade. 

Blade Bending Moment 

The blade element illustrated in Figure 1 is shown in 
Figure 2 as part of a blade section which extends from 
r, to the blade tip. Loadings due to the aerodynamic 
forces and the centrifugal accelerations are illus- 
trated in Figure 2. For a blade element of mass per 
unit span, p, rotating at a constant angular velocity, 
ft» the contribution of the above loadings on the 
element to the flatwise bending moment at position s, 
is 

dM^ = (s-s j ) [w-pv)dr - psft^ (v-v )dr ( 1 ) 

Here w is the aerodynamic load per unit span. Rigid 
body motion of the rotor blade due to teeter motion 
yields v = rsini^ = rip so that the contribution to the 
flapwise bending moment becomes 

dM ! = (r-r x )[w - U r(* + fl 2 ifO]dr (2) 

The motion of the rotor about the teeter axis when y 
Is small and the teeter axis is perpendicular to the 
blade span (63 * 0 ) is described by 

l(v+ ft 2 *) = M net (3) 
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where 


2 

1=2 / r ydr 

o 

and the net teeter moment is 

M » / rwdr - / rwdr 

net Blade 1 Blade 2 


necessary for the MOD-2 rotor to fully respond to a 
step change in wind conditions. Both the wind shear 
and the yaw produce cyclic changes in the angle of 
(4) attack. Wind shear changes the velocity component 

normal to the blades, while the yaw produces cyclic 
changes in the component of the relative velocity 
parallel to the blade. Thus, shear-induced changes in 
angle of attack are greatest at the blade tips while 
yaw cyclic angles of attack are greatest near the 
blade roots. 


Substitution of Eq. (3) into Eq. (2) yields an expres- 
sion for the differential flapwise contributions of 
the aerodynamic forces and the accelerations (both 
centrifugal and teetering) to the flapwise bending 
moment . 

H 

dM x = (r-rjwdr - J p- (r-r^rydr (5) 


Measurements of the flapwise bending moment of the 
MOD-2 wind turbine were made at station 370. At this 
station, the contributions given in Eq. (6) can be 
written as 


[aero 

Gravitational contributions must also be included. 
Including the blade twist at station 370 (twist - 0 - 
2.5°) the bending moments due to gravity are 


g / (r-rj ydr [sin^cosBsin© + sin0cos0] (6) 
r l 

For the MOD- 2 with small ^ the above contribution is 

3,219,000 iJjsinS + 142,900 cos6 N-m 

The total of the gravitational, aerodynamic, and 
acceleration contributions to the flapwise bending 
moment at station 370 is 


Aerodynamics 

The aerodynamic forces on the blades were obtained 
from integration of the local forces obtained from 
strip theory. Equating the blade forces on a blade 
element to the momentum change in differential ele- 
ment, the induced axial velocity was calculated at 
each angular position at stations along the blade 
radius. The net aerodynamic moment causing the rotor 
to teeter was calculated and used to determine the 
teeter angular velocity. Integration of the teeter 
angular velocity, yielded the teeter angle. The 
blade forces were determined using the velocity of the 
air relative to the blades so that the local teeter 
velocity, r cos 4^* (th e *®D-2 has no 60 played a 
significant role in determining the induced axial 
velocity of the air. The tangential induced velocity 
was neglected. 

Strip theory uses the concept of a fully developed 
wake in the calculation of the induced axial velocity 
so that there is an implicit assumption that an equi- 
librium-wake exists. That is, the wake is always in 
equilibrium with conditions on the blade. 

Static aerodynamic data for the NACA 23000 series 
airfoil used on the M0I>-2 were obtained from reference 
[5] . No correction was included in the calculations 
for the gap between the movable tips and the fixed 
portion of the blade and the Prandtl tip loss factor 
was used to account for decay of the circulation at 
the blade tips. 


Dynamic Stall Jfodel 


M 


1 


Total aero 


-0.352M „ + 3.219 • 10 4«in9 

net 


+ 1 .43 • 10 cos0 


(7) 


The Gormont approach [6] was selected to model the 
effects of dynamic stall based on the simplicity of 
the method and success in treating dynamic stall in 
vertical axis wind turbine applications [7], In the 
Gormont method 


The gravity terms play a large role in determining the 
cyclic bending moments. 

Wind Input 

The wind input used to model the effects of dynamic 
stall consisted of a steady wind with a yaw angle of 
15° and with wind shear. The surface roughness used 
was 0.32 m which yields a peak wind speed of 1.107 
V. , when the blade tip is at its highest position and 
a minimum wind speed of 0.736 V hub when the blade top 
is at its lowest position. Examination of Goldendale 
test site wind measurement records indicated that the 
mean shear corresponds to a surface roughness of 
0.32 m. This surface roughness yields a wind profile 
that is close to a power law profile exponent of 0.2. 

The yaw angle of 13° was selected after examination of 
the PG&E MOD- 2 wind turbine data. It was observed 
that yaw angles of 13° were sustained for periods in 
excess of 30 seconds, which is the length of time 


°L 


“r" V 


static 


( 8 ) 


where 


“r “ “static 



[.25 + 



(9) 


and Yt i s a dynamic stall parameter that depends upon 
the airfoil shape and Mach number. Since the MJD-2 
blade varies considerably in thickness ratio along the 
span, varying from 12% thickness ratio at the tips to 
27.2% thickness ratio at r * 0.3R, there was little 
data available on the parameter y l * Using extrapola- 
tion of the data given in Gormont* s report [6], the 
following form of y^ was ob tained and used for thlS 
study: 
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\ ■ ( 2 * 66 - 6 ’ 98 ^ c ) t 1 - 1.035 - 2.5 t/c J 


( 10 ) 


Table 1 below gives the variation of over the MDD-2 
blade using a local tfech number, M = rft/a. 


Table 1. Variation of with Radial Station 


r/R 

t/c 


1.00 

.12 

1.21 

.85 

.17 

.98 

.70 

.22 

.74 

.50 

.24 

.68 

.30 

.27 

.60 


The extrapolation used to characterize the dynamic 
stall parameter does not warrant the display of two- 
digit accuracy shown in Table 1; the values are shown 
to indicate the magnitudes used in this study. 

Cases Examined 

The effect of dynamic stall on the mean flapwise and 
cyclic bending moments was examined at three wind 
speeds, 9, 11.5, and 14 m/s. Cases were run with and 
without dynamic stall. The calculations were started 
with the rotor in horizontal position with zero teeter 
angle and the history was determined for six revolu- 
tions. Calculations at 14 m/s were made for eight 
revolutions. The largest cyclic loads were obtained 
at the final revolution with no appreciable 
differences being observed between the sixth and 
eighth revolution. Integration increments of 0.05 R 
and 0.10 R were used and it was found that the smaller 
integration increment resulted in an 8.7% increase in 
mean flatwise bending moment and an increase of 0.8% 
in the cyclic bending moment. 

The results reported are for the 0.10 R integration 
increment, the larger integration increment having a 
significant effect on the computation costs, approxi- 
mately cutting the costs in half. 

All the reported runs were made with fixed pitch of 2° 
on the partial span control surfaces. Interaction 
between the wind conditions and the control system was 
not Included in this study. 

RESULTS 

The results obtained for the M0D-2 wind turbine are 
shown in Figures 3 and 4. Figure 3 shows the mean 
flapwise bending moment measured at the Goldendale, WA 
test site as a function of the mean wind speed meas- 
ured at the BPA tower. Also shown In Figure 3 is the 
mean flapwise bending moment calculated in this 
study. There is good agreement between the calculated 
mean flapwise bending moment and the measured value. 
The mean flapwise bending moment is almost exclusively 
due to aerodynamic contributions since the net 
contribution from gravity and teeter Is zero. 

Figure 4 shows the cyclic flapwise bending moment data 
for Goldendale Unit #3. The data shown for 50% occur- 
rence are based on measurements, and the 99.9% case is 
based on calculations using the 50% data. 

The cyclic moments shown as the 50% occurrence were 
obtained by statistical analysis of the amplitudes of 
bending moment variations from the mean taken over 
many cycles. Another method of representing the 


cyclic variations is to take the cyclic amplitude as 
1/2 ( M max - ^j,£ n ) where is the maximum value in a 
cyclic and M tnin is the minimum value of the cycle. 
Table 2, shown below, gives the results in tabular 
form. Cyclic flapwise bendings calculated using the 
1/2 ( M tnax - *^ n ) method are given along with (ft &ax “ 
M) and the test results. Using either approach, the 
calculated cyclic bending moment is below the measured 
value . 

Dynamic stall for the cases considered is observed to 
have a small effect. At the lowest wind speed consid- 
ered, 9 m/s, there was virtually no effect of dynamic 
stall, due to the fact that angles of attack were 
below static stall over most of the blade. As the 
wind velocity was Increased, dynamic stall effects on 
the bending moment increase, but the magnitudes of 
these increases were modest. The results at 14 m/s 
show that dynamic stall increased the mean bending 
moment by 2.5% and the cyclic bending moment Increased 
by 3.3% to 9.0% depending upon the method used to 
calculate the cyclic bending moment. Regardless of 
the method used to calculate the cyclic bending mo- 
ment, the calculated values were less than the measur- 
ed values at the cyclic bending moment obtained from 
reference [3] • 


Table 2. Calculated Mean and Cyclic Flapwise 
Bending Moments at Station 370 


Wind 



Xa>T 


AM , . ** 
cyclic 

Speed 

Dynami c 

M 

W* 

W M 

(Data) 

m/s 

Stall 

N*m 

N*m 

N *ra 

N*m 

9.0 

No 

1929. 10 8 

172317 

202730 

241000 

9.0 

Yes 

1932 -lO 3 

173656 

205033 

11.5 

No 

2455*10;* 

193661 

232386 

299000 

11.5 

Yes 

248 1 *10 3 

198066 

248431 

14.0 

No 

2848 * 10;? 

229245 

274480 

357000 

14.0 

Yes 

2918* 10 3 

236758 

299075 


* Plotted in Figure 3 

**Cyclic flapwise bending moment data, 50% occurrence 
from Figure 4 


CONCLUSIONS 

The results obtained for cyclic and mean flapwise 
bending moment for cases with and without dynamic 
stall show no substantial change due to the Inclusion 
of dynamic stall. Further comparisons of calculated 
and measured bending moments show good agreement for 
the mean bending moments and poor agreement for the 
cyclic bending moments. The cyclic bending moments 
calculated were for a 15° yaw angle which is a 
frequent, but not the usual, case for the M0D-2. It 
is concluded that spatial variations in wind speed due 
to wind shear and yaw are not the primary drivers of 
cyclic loading. Further, the effect of dynamic stall 
on the mean and cyclic bending moments is small. Wind 
turbulence and rotor flexibility were omitted in this 
preliminary study of the effects of dynamic stall and, 
of these omissions, wind turbulence is felt to play a 
significant role in determining the magnitude of the 
cyclic loads. The role of dynamic stall in interact- 
ing with wind turbulence remains unknown with the 
major barrier to understanding being the inclusion of 
wind turbulence in the loading analysis. 
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The following Appendices involve subjects of general 
interest to HAWT developers. 

APPENDIX A 

WAKE EXPANSION MODEL 

A simple analytical model for the wake and flow field 
of a horizontal axis wind turbine can be developed 
from actuator disk theory. Let_ the mean axial velo- 
city at the rotor by u * V(l-a) and tjie mean axial 
velocity in the far wake by u fW - V(l-Xa). The axial 
velocity along the axis of a rotor is approximated by 


which yields a linear variation with r for the radial 
velocity, v . The linear variation in radial velocity 
is in qualitative agreement with the results of a free 
wake analysis. 


Streamlines in 

the flow 

may be obtained 

f rom the 

relation: 




r— ) 

(— ) 

v 

r 

(A-3) 

Mx ' 

0 S.L . 

^ S.L. 

v 

X 



Integration of this equation yields 


Rwake * fi - fi f 1/2 (A-4) 

R 1+q ' 


where E = a(X-l)/(l-a) . Equation (A-4) may also be 
used for any streamline that starts at the rotor. 

In order to use the above approximation relation, 
values of a and X are required. The value of a is a 
mean value of the induction at the rotor and can be 
approximated from the thrust coefficient. 

For C T < 0.64, X = 2 and 

* 4a (1-a) (A-5) 

which is a quadratic equation. 

For C^ > 0.64 

C T > 0.16 + 2.4l 0-6) 

while X can be obtained from 



2a( 1-a) 


The equation for the wake boundary is an approximation 
and although it satisfies continuity, Eq. (A-4) has 
positive curvature for the streamlines at the rotor. 
The point of zero curvature for the streamlines is 
defined by 


(X-l)q 

v x = v[l-a - a / 1+q 2 ] 


O-l) 


where q - x/R, R is the rotor radius and the rotor is 
at the origin with x measured downstream. The above 
relation is exact for the velocity along the axis of a 
rotor with a cylindrical vortex sheet wake when X 
2. Using the above equation for v x throughout the 
flow field one obtains from the continuity equation 


v 

r 


Va( x-l)n 

2(.*<, 2 ) V2 ' 


r 



0 - 2 ) 


E(q 2 + 5/6) = q /1+q 2 


The point of zero streamline curvature varies from the 
rotor at E=0 (no wake expansion) to q ** 1 .42 for Ef=l 
(infinite wake diameter) . 

APPENDIX B 

FREE YAW OF RIGID HUB ROTORS 

The stability of rigid hub rotors in free yaw may be 
described by examination of the change in yaw moment 
with the yaw angle, 6. 

3M 

k = (B-l) 

* " 35 
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Other yaw moments can cause yaw tracking errors but 
the yaw stability is determined by magnitude and the 
algebraic sign of k. While yaw tracking errors are 
mostly due to the imbalance of in-plane aerodynamic 
forces due to wind shear and tower shadow, the spring 
constant, k, is due to the normal aerodynamic for- 
ces. Figure B-l illustrates the role of the coning 
angle ^ in the stability of rigid hub rotors in free 
yaw. 

A rotor with blade horizontal is shown schematically 
when viewed from above. The velocity normal to the 
blades is shown for blade elements on the right hand 
side and the left hand side of the rotor. The veloci- 
ty normal to the blade determines the angle of attack; 
that is 

a ~ (B-2) 

The induction may be different on opposite sides of 
the rotor; however, the primary dependence is due to 
the normal velocity. It may be observed that 


\ > \ 

hence the angles of attack for the illustrated blade 
will obey 


°R > “l 

for the illustrated figure. 


CBr-3) 


The normal forces on the right hand and left hand 
sides are proportional to a Cl , hence the moment 
about point 0 tending to yaw ^he rotor in the <$- 
direction is proportional to 

-Cl (ag - «l) Cfr-4) 
a 

The spring constant, k, is the result of integration 
of the differential yaw moment due to yaw over the 
rotor. If Cl >0 the illustrated rotor will be 
stable In yaw a since the coning angle ij; produces a 
greater normal force on the blade occupying the right 
hand side of the rotor than on the left. However, 
beyond stall for many airfoil sections, C La < 0 , and 
the rotor yaw moment will Increase with increasing 
yaw. In such a case the rotor will operate In stable 
mode at an upwind position where the coning angle is 
negative. Thus the stability of rigid rotors horizon- 
tal axis wind turbines In yaw can be characterized by 
the integrated contribution of Cl over the rotor* 
For a rotor with three or more blades (B = number of 
blades) the dominant term in the yaw "spring constant" 
is 


1 ? 

; (f) n x Cl 4 -cos *dn 

o a 


where 


n = r/R 



<j> is the local flow angle with respect to the 
rotor disk. 



Figure 1. Rotor Geometry and Coordinates 



Figure 2. Blade Section with Aerodynamic 
and Centrifugal Loading 
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ROTOR 


Figure 3. Mean Flapwise Bending Moment at Station 
370. Data from Reference [3]. 


Figure B-l. Wind Turbine in Yaw. 



Figure 4. Cyclic Flapwise Bending Moment at Station 

370. Data from Reference [3]. Calculated 

Cyclic Bending Moment is ~ M ). 

J max min 
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